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Interrelation Between Safety Factors and Reliability: Part 1 
Random Actual Stress & Deterministic Yield Stress 


Isaac Elishakoff 

Department of Mechanical Engineering 
Florida Atlantic University 
Boca Raton, FL33431 


'The factor of safety was a useful invention of the engineer a long time 
ago that served him well. But it now quite outlived its usefulness and 
lias become a serious threat to real progress in design” 

D. Faulkner 

“Current structural safety design practices are considered inadequate 
for future launch vehicles and spacecraft.” 

V. Verderaime 

“[Deterministic safety measures] ignore much information which may 
be available about uncertainties in structural strengths or applied 
loads.” 

R. Melchers 


1. Introduction 

Attempts at probabilistic interpretation of the deterministic safety factor have 
been made in the literature despite the fact that the “spirits” of these two approaches are 
entirely different. Before discussing them, it is instructive to quote some representative 
excerpts from popular textbooks concerning its definition: 

(a) “To allow for accidental overloading of the structure, as well as for possible 
inaccuracies in the construction and possible unknown variables in the analysis of 
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the structure, a factor of safety is normally provided by choosing an allowable 
stress (or working stress) below the proportional limit”. 

(b) “Although not commonly used, perhaps a better term for this ratio is factor of 
ignorance”. 

(c) “The need for the safety margin is apparent for many reasons: stress itself is 
seldom uniform; materials lack the homogeneous properties theoretically assigned 
to the abnormal loads might occur; manufacturing processes often impart 
dangerous stresses within the component. These and other factors make it 
necessary to select working stresses substantially below those known to cause 
failure”. 

(d) “A factor of safety is used in the design of structures to allow for (1) uncertainty 
of loading. (2) the statistical variation of material strengths, (3) inaccuracies in 
geometry and theory, and (4) the grave consequences of failure of some 
structures”. 

(e) “Factor of safety (Nfs), where Nfs > 1 , is the ratio of material strength (usually 
ultimate strength or yield point) to actual or calculated stress. Alternatively, factor 
of safety can be defined as the ratio of load at failure to actual or calculated load. 
The factor of safety provides a margin of safety to account for uncetainties such 
as errors in predicting loading of a part, variations in material properties, and 
differences between the ideal model and actual material behavior” (Wilson, 
1997). 

(f) “Choosing the safety factor is often a confusing proposition for the beginning 
designer. The safety factor can be thought of as a measure of the designer’s 
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uncertainty in the analytical models, failure theories, and material property data 
used, and should be chosen accordingly . . . Nothing is absolute in engineering any 
more than in any other endeavor. The strength of materials may vary from sample 
to sample. The actual size of different examples of the “same” parts made in 
quantity may vary due to manufacturing tolerance. As a result, we should take the 
statistical distributuions into account in our calculations.” (Norton, 2000). 
Freudenthal remarks [20] “ ... it seems absurd to strive for more and more 
refinement of methods of stress-analysis if in order to determine the dimensions of the 
structural elements, its results are subsequently compared with so called working stress, 
derived in a rather crude manner by dividing the values of somewhat dubious material 
parameters obtained in conventional materials tests by still more dubious empirical 
numbers called safety factors”. 

Indeed, it appears to the present writer that in addition to its role as a “safety” 
parameter for the structure, it is intended as “personal insurance” factor of sorts for the 
design companies. 

2. Four Different Probabilistic Definitions of a Safety Factor 

Consider an element subjected to a stress a. Let it be a random variable, denoted 
by capital letter Z, whereas a lower case notation describes the possible values crthat the 
random variable I may take. The strength characteristics say the yield stress may also be 
designated by upper case notation with a y being the possible values h y may take. The 
various possible definitions of the safety factor 5 are 
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( 1 ) 


£(£ v ) 

s. = — — 

E(l) 

which is referred to as the central safety; E{L y ) denotes the mathematical expectation of 
while E(L) is associated with the mathematical expectation of L 
On the other hand one can treat the ratio 


0 = ^ 
~ I 


as a random variable. Its mathematical expectation 


( 2 ) 


s 2 = E(Q) = B 


fl. 




( 3 ) 


could be also interpreted as a safety factor. The third possible definition of the safety 
factor is 


5 3 = E( l y )E 


rn 




( 4 ) 


In specific cases some of the above safety factors may coincide. For example, 
when the yield stress is random, but the stress is a deterministic quantity, i.e. takes a 
single value a with unity probability, we have 


5 , = 5 , = 5 , 


If stress is random, but yield stress is a deterministic quantity a y , then 


5, = S, 


( 5 ) 


( 6 ) 


The fourth definition of the safety factor was proposed by Birger (1970). He 
considered the probability distribution function of the random variable Q 
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( 7 ) 


f% "| 

F Q (q) = Prob -r < q 

V L J 

Then he demands this function to equal some value po. 

(L \ 

Fg(q) ~ Prob -=r^q = p 0 (8) 

L y 

The value of q = qo that corresponds to the /-> 0 th ffactile of the distribution function Fg(q) 
is declared as the safety factor. This implies, that of po = 0.01, and say qo ~ 1, that in 
about 99% of the realizations of the structure the deterministic safety factor will be not 
less than 1 3 . 

We are asking ourselves the following question: Can we express the safety factors 
by probabilistic characterization of the structural performance? The central idea of the 
probabilistic design of structures is reliability, i.e. probability that the structure will 
perform its mission adequately, as required. In our context the mission itself is defined 
deterministically, namely we are interested in the event 

2 < 2 , ( 9 ) 

i.e. that the actual stress is less than the yield stress. 

Since both Z and £ v may take values from a finite or infinite range of values, the 
inequality (9) will not always take place. For some realizations of random variables Z and 
Z y the inequality may be satisfied, whereas for the other ones it may be violated. 
Engineers are interested in the probability that the inequality (9) will hold. Such a 
probability is called reliability, denoted by R 

R = Prob(Z < Z y ) (10) 

Its complement 
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P f = l-f? = Prob(£>E y ) 


( 11 ) 


is called the probability of failure. It is a probability that the stress will be equal to or will 
exceed the yield stress. It is understandable that engineers want to achieve a very high 
reliability, allowing, if at all, a extremely small probability of failure. 

It appears, at the first glance, that the approaches, based on the deterministic 
allocation of the safety factors, or that based on reliability design are totally 
contradictory. We will pursue this subject in more detail. In this report and its companion 
(report #2), we discuss particular cases, whereas at a later stage (in our report #3) we will 
pursue the general case, in which both I and £,< will be treated as random variables. 


3. Case 1: Stress Has an Uniform Probability Density, Strength 
Is Deterministic 

Let the stress £ be a random variable with the uniform probability density 


/i(*H 


-, for a L < a <a u 


( 12 ) 


0, otherwise 


where <jl is the lowest value that the stress may take, whereas <ju is the greatest value the 
stress may assume. We treat the yield stress £>. to be a deterministic quantity, i.e. to take a 
single value Oy with unity probability. The probability distribution function 


U 

F z {a) = Prob(I < ct) = J ffafda 


(13) 


reads 
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0, for cr < a. 


F z (g) = \ 



for <j l < a < g v 


[ 1 , for G a < G 


(14) 


The reliability reads 


R = Prob(E < I. y ) - Prob(Z < g } .) (15) 

or, in light of Eq. (13) we get 


R = F t (G y ) (16) 

In other words, the reliability equals the stress distribution function F z evaluated at the 
yield stress (Fig. 1). 
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Fig. 1 Reliability equals the probability distribution function of the 
actual stress evaluated at the level of yield stress 
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Bearing in mind Eq. (14) we get 


[0, for <J y <o L 


R = 


a v - cj l 

< — , for <j L < o y < o v 

a L , a L 

}, for cr L , < o y 


( 17 ) 


This formula can be rewritten in more convenient form. We note that the mean value of 
the stress equals 

£(!) = \{o L +o v ) (18) 


whereas the variance of the stress equals 

Var(Z) = ; 2 -(o v - o L ) 2 


(19) 


From these two equations we first express the denominator in Eq. (17) 


cry - o L - yj\2Var(I,) 

as well as the lowest possible value the stress can take 

2£(L) - Vl2Far(Z) 

<y j — 

L 2 
= E(L)-y[ 3Far(I) 


Let 


( 20 ) 


( 21 ) 


o L <o <o u (22) 

then, in accordance with Eq. (17), we have 

a -£(!)- J 3Ehrr(X) 

R= _y w ( 23 ) 

2^3Var(l) 

By dividing both the numerator and the denominator by the mean stress E(L), and 
introducing the coefficient of variation of the actual stress 
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(24) 


v i = 


VFar(S) 

E&) 


we get, instead of Eq. (23) 


5, — 1 — V3v £ 

2a/3v t 


(25) 


As is seen reliability is directly expressed in terms of the central safety factor S\ and the 
coefficient of variation of the involved random variable vj. Thus, the reliability methods 
allow to rigorously , rather than arbitrarily introduce the safety factors. The safety factor 
is expressed from Eq. (25) as follows 

5, =1 + V3v s (1 + 2/?) (26) 

Maximum value of the safety factor is achieved when the reliability tends to unity from 
below 


< 27 > 

For example, for coefficient of variation 0.05 the safety factor assumes the value 1.26; for 
the coefficient of variation 0.1 the safety factor equals 1.52; for coefficient of variation 
0.15 it takes a value 1.78 etc. We conclude that with greater variation of the involved 
random variable, the safety factor must be increased This qualitative conclusion is in 
line with our anticipation. Yet, it is seen that the reliability context allows one to make 
quantitative judgements in terms of the required reliability and the coefficient of 
variability. 
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4. Case 2: Stress Has an Exponential Probability Density, Yield 
Stress Is Deterministic 


Consider now that the stress has an exponential probability density 

[0, for a< 0 

A(°) = i 

[aexp(-aa), for a > 0 
The corresponding probability distribution function reads 

Fv(cr) - Prob(S < cr v ) = [1 - exp(-acr)](/(cr) 


(28) 


(29) 


where U(o ) is the unit step function, it equals unity for positive y and vanishes otherwise. 
The parameter a is reciprocal to mean value of stress 

E(l) = - (30) 

a 

Also, since parameter a is the only free parameter in the density (28) all probabilistic 
moments depend solely upon it. Thus, variance also is expressible in terms of a , as 
follows: 


Var(Z) 


(31) 


a 


Since the coefficient of variation 


Vfar(I) _ 1 !_a 
Vz ~ E(Z) ~ M a 


(32) 


is unity, or 100%, we must anticipate high levels of safety factor, in order to ensure the 
high level of required reliability. The latter equals, in view of Eq. (19) 

R - Pro6(I < cr v ) 

= [1 - exp( -a cr y )]£/ (cr y ) 

We first express a from Eq. (30) as 


( 33 ) 
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a = !/£(!) 


(34) 


and substitute it into Eq. (33), to arrive at 

V(a y ) (35) 

In view of the central safety factor si, Eq. (35) is rewritten as 

/? = [l-exp(-5)]f/(cr > .) (36) 

As is seen, a direct relationship is being established between the safety factor and the 
reliability. Once the required reliability is specified the associated safety factor equals 

s = In — - — (37) 

1 - R 

For example, reliability of 0.9 leads to the safety factor 2.3; the reliability of 0.95 results 
in safety 3 etc. Such high values, as indicated above stem from the fact that the stress 
exponential probability density is associated with high, namely 100% variability. It is 
immediately seen, that one of the reasons for the high variations in this particular case is 
the fact that the stress can take any value on the positive axis. 


R = 


1 - exp 


£(S) 


5. Case 3: Stress Has a Rayleigh Probability Density, Yield 
Stress is Deterministic 


Consider now the case in which the stress is characterized by a Rayleigh 
probability density: 




f0, for a < 0 

for o> 0 


(38) 
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The approximate probability distribution function is 


u 

F z (o) = Prob(E < <r) = J f z (a)da = 


1 - exp 


f 2 \ 

' a ' 


V 2 


\U(a) 


(39) 


The reliability evaluation reads: 

R = Prob(I < u y ) - F z (cr y ) (40) 

Again, the reliability equals the probability distribution function of the stress evaluated at 
the level of the yield stress. Hence, in view of Eq. (39) 


R = 


1 - exp 


v 21,2 J 


U(CT y ) 


(41) 


We would like now to express parameter b in Eqs. (38) and (39) through the probabilistic 
characterization of the stress: 


oo r J 

E( I) = jqf z (o)da = * 1.256 


(42) 


Var(l)= ](a-E(l)) 2 Ma)da 


— b 2 * 0.436 2 


(43) 


We express b from Eq. (42) 


b * = 0.8£(L) 

1.25 


and substitute it into Eq. (41) to yield 


r 

f 

2 \] 

/T 

< 1 - exp 

— 

V 2 

1 


2[0.8£’(E)] 2 J 

f 

f 

0.781250-/^ 

« 1 - exp 

K 

[£(E)] 2 y 


(44) 


( 45 ) 
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We take into account the definition of the central safety factor to get 

R=[ 1 - exp(-O78125j/)]C/(0V) ( 46 ) 

This formula allows to express the safety factor by the reliability 

s*l. 13 Jin— (47) 

V 1 -R 

Thus, the reliability R = 0.9 yields in central safety factor 1.71, the reliability of 0.95 
results in safety factor 1.96. The required reliability of 0.99 is associated with safety 
factor 2.42 etc. Again, reason for these values is the high coefficient of variation. Indeed, 
Eqs. (42) and (43) suggest that the coefficient of variation equals: 

(48) 

£(I) 1.256 

Although this is a smaller variability than in the case of the stress with exponential 
probability density, still, hopefully, 52% variation is seldom encountered in practice. 


6. Case 4: Stress Has a Normal Probability Density, Yield Stress 
Is Deterministic 

We consider now the case in which the stress is characterized by a normal 
probability density 

/l<CT) = wt exp 

The distribution Sanction reads 


1 f cr -a 

2 


00 < or < 00 


(49) 


J 
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1 ^ 

/v(o-) = —7== f exp 
b\2n _ J „ 


1 

2l * ; 


k* 


= O 


^ <7 -a'' 


1 ^ 


k* 




where a is the mean stress, and £ is the mean square deviation, 

E(Z) = a 
Var(L) = b 2 


The reliability equals 

( a - a ^ 

i? = Pr ob(S < <7 V ) = ) = O 1 


(50b) 


(51) 

(52) 


(53) 


or with Eq. (50) taken into account 


R = <& 


Vy-EV) 


(54) 


JVar(I.) 

Dividing both the numerator and the denominator by £(L) we rewrite Eq. (53) as follows: 

R = Ol 2_l | (55) 


K v z J 

where si is the central safety factor, vv is the coefficient variation of the stress. Eq. (54) 
allows again to express the safety factor via the reliability 

5] = l + v s O'(7?) (56) 

where 0 _1 (7?) is a function that is inverse to 0(7?) 

We note the following values of the inverse normal probability function (Benjamin and 
Cornell, 1970, p. 655) 
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0“' (0.9) = O' (1-10“') = 1.28 
O'(0.99) = 0“'(1 -10“ 2 ) = 2.32 
0“' (0.999) = 0“’(1 - 10“ 3 ) = 3.09 
O“'(0.9 999) = 0“’(1 - 10 4 ) = 3.72 
O '(0.99 999) = O'(l-10‘ 5 ) = 4.27 

O' (0.999 999) = 0> '(1 - 10“ 6 ) = 4.75 (57) 

$ '(0.9 999 999) = O' (1 - 10“ 7 ) = 5.20 
O ‘(0.99 999 999) = O“‘(l-10“ 8 ) = 5.61 
O ‘(0.999 999 999) = 0"‘(1 - 10“ 9 ) = 6.00 
O“‘(0.9 999 999 999) = 0“'(1 - 10 10 ) = 6.36 
O '(0.99 999 999 999) = O ‘(1 - 10“") = 6.71 

Thus, the safety factor becomes, for the coefficient of variation equal 0.05, respectively 

s, =0.050 ’(0.9)= 1.064, for 7? = 1-10‘ 

5, =0.050 ’ (0.99) = 1.116, for tf = l-10“ 2 
5, =0.050 ■'(0.999) = 1.155, for 7? = 1-10“ 3 
5 , = 0.050“' (0.9 999) = 1 1 86, for R = 1 - 10“ 4 
5 , = 0.050 ' (0.99 999) = 1.2135, for 7? = 1 - 1 0“ 5 

5, =0.050-' (0.999 999) = 1.2375, for 7? = 1-10“ 6 (58) 

5 , =0.050 ‘(0.9 999 999) = 1.26, for 7? = 1-10 7 

5 , =0.050“' (0.99 999 999) = 1.2805, for fl = l-10“ 8 

5 , =0.050“' (0.999 999 999) = 1.3, for 7? = 1-10“ 9 

5 1 = 0.050 ‘(0.9 999 999 999)= 1.318, for R = 1 — 10 10 

5 , =0.050“' (0.99 999 999 999) = 1.3355, for R = 1-10“" 

As is een, there is direct relationship between the safety factor and required reliability. 
One can suggest asymptotic relationship between the safety factor and reliability. One 
observes from Eq. (57) that the knowledge of the coefficient of variation (the ratio 
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between the standard deviation and the mean) and required reliability directly yields the 
level of the required safety factor. 


7. Case 5: Actual Stress Has a Log-Normal Probability Density, 
Yield Stress is Deterministic 

Consider now the case in which the actual stress E is distributed log-normally, with 
the following probability density function: 


M<y) = 


1 


exp 


(In a -a t )' 


2b y 


a > 0 


o b, -Jin 

and vanishes otherwise. The mean value of the stress equals 

E(L) = exp(a z + \b^) 

whereas the variance reads: 

Var(L ) = exp(2a s + 6 I 2 )[exp(Z> £ :! )-l] 

The reliability equals 

R = Prob(I < cr y ) = Fv(cr v ) 

The probability distribution function for the log-normal variable E is 




We make a substitution 


f 1 

1 pyh 



2 



x j 



dt 


In t - a T 


( 59 ) 


( 60 ) 


( 61 ) 


( 62 ) 


( 63 ) 


= z 


( 64 ) 


to obtain 
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( 65 ) 


b Z 1 

F z(&)= | -7==exp(-ir^ 


In <j - a v 


In view of Eq. (65), expression for the reliability 


R = <& 


In cr v - a z 


Central safety factor equals 




E(l) exp (a z + -b z ) 


Knowing parameters a z and b z determines both the central safety factor and the 
reliability. 

On the other hand, if E(L) and Var(L ) are given, one needs the formula of 


transformation from Eqs. (60) and (61). We substitute Eq. (67) into Eq. (66) to get 


R = <& 


In <7 V - In E( I) - - }n[£ 2 (E) + Var(l)] - ln[£ 2 (I)]f 


/ln[£ 2 (E) + Var(Z)] - ln[£ 2 (I)] 


Thus, the find formula can be rewritten as: 


„ Jln5,-i[ln(lW)] 

I\ = ===== 

yjlnil + Uz 2 ) 


We are unaware of the other derivation of this expression. 
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8. Case 6: Actual Stress Has a Weibuil Probability Density, 


Strength is Deterministic 


Let us study the case of the probability of the distribution function stress 

Fv{<j) = exp - exp 


r Or - <7^ 


\ j 


(70) 


The reliability, therefore, is given by 

R = F z( a v) = ex P 


exp 


f — o y ^ 


V u i J J 


(71) 


According to Haidar and Mahadevan (2000) who do not deal with the material in this 
section, but use the Weibuil distribution, the mean value and the variance can be 
expressed via a z and b z analytically. In our setting their formulas read: 

1 1 


K 


Thus, 


b z Vct yjVaril) 
ct z = £(!)- 0.5772k 

E(L) = ci z +0.5772*. 

2 

Var(L) = —b 2 
6 


(72) 


(73) 

(74) 


Reliability becomes: 


R = exp 


exp 


' E(Z)-o.sn2^^- a , 

71 ' 


V6For(I) 


n 


J J 


(75) 
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or, dividing both numerator and denominator by E(L) and recalling definition of the 
central safety factor $, = a y / E(L) and of the coefficient of variation of the actual stress 


v z = y JVar('L)/E(I) we get 

R = exp 


exp 


^1 -0.45v e -$, ^ 
y 0.78v s j 


(76) 


This formula too apparently is given for the first time. It connects the reliability with the 
central safety factor $i and the variability v z . Conversely, if the required reliability is 
specified, one can directly determine the safety factor 

7 

(77) 


$, = 1 - v r 

0.78 In 

r, n 

-0.45 

i L 


l Rj 



The following values are obtained for v E = 0.05: 

R = 0.9, $,=1.11 

R = 0.95, $, =1.14 

R = 0.99, $,=1.18 

/? = 0.999, $,=1.29 

R = 0.9 999, $, =1.36 
R = 0.99 999, $, = 1.45 

For v z = 0.1 we get 

R = 0.9, $,=1.22 

R = 0.95, $,=1.28 

R = 0.99, $,=1.36 

R = 0.999, $,=1.58 

R = 0.9 999, $,=1.72 
R - 0.99 999, $, =1.90 


(78) 


(79) 
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9. Conclusion 


As is observed from this report the use of the safety factor is not contradictory to 
the employment of the probabilistic methods. Moreover, in many cases the safety 
factors can be directly expressed by the required reliability levels. However, there is a 
major difference that must be emphasized: whereas the safety factors are allocated in 
an ad hoc manner, the probabilistic approach offers a unified mathematical 
framework. The establishment of the interrelation between the concepts opens an 
avenue for rational of safety factors, based on reliability. 

If there are several forms of failure then the allocation of safety factors should be 
based on having the same reliability associated with each failure modes. This 
immediately suggests, that by the probabilistic methods the existing overdesign or 
underdesign can be eliminated. 

This is done by calibration of the reliability levels with one of the safety factors 
that is already accepted. Thus, via such an approach, the other failure modes’ safety 
factors can be established 

This is illustrated fig. 2, which shows that presently safety factor are assigned in 
an ad hoc manner to each failure mode, but there is no interrelation between them. 
Fig. 3 illustrates the consistent allocation of the safety factors can be performed. 

The report No. 2 deals with the reverse case, namely when the actual stress is 
deterministic, but the yield stress is random. Report No. 3 will discuss the general 
case in which both the actual stress and the yield stress are treated as random 
quantities with the attendant interrelationship between the reliability and safety 
factors 
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Safety factor n ( ^ 



Fig. 2 Present Status: No Connection Between Safety Factors, Leading 

to Overdesign or Underdesign 
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Fig. 3 Future Status: Equal Reliability Allocation May Connect Dafety 

Factors 
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Interrelation Between Safety Factors and Reliability: Part 2 


Deterministic Actual Stress & Random Yield Stress 

I Elishakoff 

Department of Mechanical Engineering 
Florida Atlantic University 
Boca Raton, FL 33432 

In the previous report we studied the case in which the actual stress was treated as 
a random variable, while the yield stress was considered as a deterministic quantity. In 
this report we investigate the reverse case, namely, when the actual stress is deterministic, 
while the yield stress is treated as a random variable. Various probability densities to 
model the actual behavior of the structural element in question are considered. 

Case 1: Yields Stress Has an Uniform Probability Density, Actual 
Stress is Deterministic 

Let the yield stress have an uniform probability density 

f 1 

M cr y L < a < a v U 

A y i. a y)-\°y- v ° vL ( 1 ) 

0, otherwise 

L 7 

where a vL is the lower possible level the yield stress may take; a y V is the upper 

possible level the yield stress may assume. 

The probability distribution function of the yield stress reads 
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0, for o y < o yL 




°y G v,L 


for o yX <o y < o y U 


r y,U 

for cr y U < <7 y 


( 2 ) 


The reliability equals 


R - Prob(I < £ v ) = Prob(<7 < 

= Prob(S y > <t) = 1 -Prob(E v < a) 


(3) 


Thus, in view Eq. (2), we get 


R = 


1, for 

° ~ °y,l 


°y < CT v.I 


1- 

O y,U 

0, foi 


~°y,L 


for 


°y,L ^Oy<<*y,V 


(4) 


>r a y U < a y 


Consider the case in which the yield stress belongs to the interval \oy,L,o y ,v\. In this case 
from Eq. (4) we have for the reliability 


R 


< fyV_ 

~rr yL 


(5) 


We note that the mean yield stress equals 


£(£ y ) = +Vy,u) 


( 6 ) 


whereas the variance of the yield stress reads 

Var(l y ) = ^(a yU -a y _ L ) 2 (7) 

We express upper level of the yield stress Oy,u as follows, in terms of the mean yield 
stress E(L y ) and variance of the yield stress Var(L y ) via Eqs. (6) and (7): 
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( 8 ) 




2E(Z y ) + Jl2Var(?j 


= E( l y ) + y l3Var& y ) 


The denominator in Eq. (5) is directly expressible by the variance as 2f3Var(L y ) in Eq. 


(7). Thus, the reliability in Eq. (5) can be rewritten as 


R 


E(L y ) + fiVar(L y )-<T 
2j3Var(L y ) 


We divide both the numerator and denominator by crand express the ration 


fiVar&y) _ ^3 Var(X y ) £(£,) 
a E( S v ) a 

= V3Vv 5, 


where 


\ = JVar(l y )/E(l y ) 

is the coefficient of variation of the yield stress, 

s,= E^/o- 


is the central safety factor: 


R = 


s,(l + V3Vj. )-l 

2-s/3vj. 5, 


Eq. (13) allows to express the central safety factor as a function of the reliability: 

1 


5, = 


1 1+V3v £ (1-27?) 


(9) 


( 10 ) 


( 11 ) 


( 12 ) 


(13) 


(14) 


This equation is remarkable for the required reliability R is directly connected with the 
central safety factor s\. Thus, if the required reliability 0.9 is set, at the coefficient of 
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variation of the yield stress v Iv = 0.05 we get the level of safety factor equal 1.07; for R = 

0.99 we get si = 1.09, reliability level 0.999 corresponds to .Vi = 1.095. At the greater 
coefficient of variation, namely, that comprising 10% we get 

5, =1.16, for R = 0.9 

5, =1.20, for R = 0.99 (15) 

5, =1.21, for R = 0.999 

When the variability constitutes 20%, we obtain 

5, =1.38, for R = 0.9 

5, =1.51, for R = 0.99 (16) 

s, = 1.53, for R = 0.999 

etc. yielding greater needed safety factors with greater variability, if the demanded 
reliability level is fixed 


Case 2: Yield Stress Has an Exponential Probability Density, 
Actual Stress Is Deterministic 


Consider now the case in which Y, y is variable with exponential probability density 
but I is deterministic. Hence I takes only a single value a with unity probability. 

The probability density of reads: 


A (<%) = 


0, for <J y < 0 

aexp(-CH7 y ), for o y > 0 


Here f z ^ is the probability density of the yield stress. 

The probability distribution function of is defined as (Fig. 1) 


(17) 


Fv (a v ) = Prob(I <a y ) 


(18) 
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i.e. as a probability that E r will take values that are not in excess of any pre-selected value 


<*y- 


According to the definition of the probability distribution 


Jaw* 


(19) 


we get 




0, /or a y < 0 
1 - exp(-a<r v ). 


for 


dy > 0 


The parameter a is the reciprocal of the mathematical expectation 

£(£,) = - 

a 


The reliability reads 

R = Prob(Z < E v ) = Prob(<r < S v ) 

= Prob(E v > cr) = 1 -Prob(Z v < cr) 


( 20 ) 


( 21 ) 


( 22 ) 


In the right side of the equation (20) we recognize that the quantity Prob(Z v < cr) coincides 
with Eq. (16) when instead of o y in Eq. (16) we substitute cr. In other words Prob(Z y <a) 
equals the probability distribution function of the yield stress evaluated at the level of the 
actual stress (Fig. 2): 

Pr ob(I v < cr) = F ly (cr) (23) 

Thus, bearing in mind Eq. (18) we get 

fO, for a < 0 

Prob(E v < a) = J (24) 

[1 - exp(-ocr), for cr >0 
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Hence, the reliability in Eq. (20) becomes 


R = 1 - Prob(Z v < cr) = 1 - (cr ) 


(25) 


= 1 - [1 - exp (~a<j)]U (<x) 

where U(d) is the unit step function, i.e. U(o ) equals unity for positive cr and vanishes 
otherwise. Taking into account the relationship (21) we can rewrite Eq. (20) in the 
following manner 


R = < exp 


E& v ) 


mo) 


(26) 


We recognize the argument in Eq. (26) to be reciprocal of the safety factor. Due to Eq. 
(5), of the report 1, three safety factors coincide in this case. Hence we denote them by a 
single notation s. Thus we get the following relationship: 

R = exp(- 1 / s)U (cr) (27) 

As we see, reliability is intimately connected with the safety factor in the case under 
consideration. In fact, the safety factor can be expressed directly from Eq. (22) for a > 0. 


'“'**> V* (28) 

R 

In this particular case if the required reliability equals 0.9 the safety factor 1/7/2(0 9) is 
greater than 9! 

The results in this case, although may seem to be very surprising, are quite 
understandable. The variance of the yield stress 



(29) 


The coefficient of variation in this case 


N AS A/CR— 200 1-211 309 


35 



( 30 ) 


Jr*®,) l/o , 

c.o.v. = — — = = 1 

E(L y ) 1 la 

equals unity; i.e. there is a large variation around the mean value of the yield stress; 
hence, large safety factors are needed to achieve the required reliability levels. 

Case 3: Strength Has a Rayleigh Probability Density, Actual 
Stress Is Deterministic 

The probability density of the strength is given by 


/i„( CT y) = 1cr 


0, for o < 0 


< a 2 ^ 


rr ex P| 


v 262 , 


with parameter b 2 . The distribution function is 


, for a > 0 


F z<°y) = 


1-exp 


t £7 2 ^ 


U(aJ 


We also note that the mean strength is (Ref 1, p75) 

£(2,) = * 1 25i 
whereas the variance of the strength is (Ref. 1, p75) 

Var(Z : ) = ^ " * 0.43 b 2 


(31) 


(32) 


(33) 


(34) 


Reliability is given by Eq. (3) 

R = 1 - Pr ob(L y < or) = 1 - F Zy (<r) 

Bearing in mind Eq. (27) we get: 


( 35 ) 
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(36) 


r ^ \ 


R = 


exp 


G 

2b 5 


TOO 


yj 


Now, taking into account Eq. (33) we can substitute instead of b , 


£(I V ) 

^-1^ = 0.8£(I V ) 
1.25 ' 


(37) 


to get 


/? = f/(cr)exp 


2[0.8£(I V )] 2 J 


(38) 


or, in terms of the safety factor 


s = 


£(2 V ) 


(39) 


we obtain 


= f/(<r)exp 


0.78125 


(40) 


Safety factor s can be expressed from the reliability 


0.78125 0.8839 


s = 


(41) 


ln(l / R) Vln(l JR) 

Let the reliability be set at R = 0.99. Eq. (40) yields safety factor 8.82. This result is again 
understandable since the coefficient of variation in this case too is quite large: 


^jVard.J ~ Vo.43* 2 

E&r) * 1.25ft 


(42) 
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Case 4: Various Factors of Safety in Buckling 


It is best to start with an engineering example, first in the deterministic setting. 
Consider an element that is simply supported at its ends. It is subjected to the 
compressive load P at the ends, as well as the concentrate bending moment. M. The 
section modulus of the cross section is denoted by S. Material’s proportionality limit a p , 

as well as the yield stress cr are given. We are interested in determining the safety factor 
in 3 different regimes. 


(a) During Use of the Element Both M and P Increase Simultaneously 

In this case we have 


M P 

_ max. + _ 

max o A 

S A 


where M max the maximum bending moment 


M = 


M 


max kL 

COS — 


where 


(43) 


(44) 


* = (45) 

Me i 

A = cross sectional area. 

Since the relationship between the stresses and the load P is nonlinear, the safety 
factor risF is determined as follows. We multiply the load by »sf so as to achieve a level 
of stress equal to the yield stress. Thus, the deterministic safety factor is derived from 
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where 


k 


V El 


(47) 


Consider now the probabilistic setting of the problem. Let c r v be a random variable L v 

Then the central safety factor J'i is determined from the equation, in the manner, 
analogous to the deterministic setting, except that o y is replaced by E(L } ), and is 
replaced by si . 


E(L y ) 


S,M 


kL 
S cos^ — 
2 


S \P 

~A 


(48) 


where 


k 


y 



(49) 


(b) During Use of the Element the Axial Force Remains Constant, 
Concentrated Moment Increases 

Deterministic safety factor hsf is found from the equation: 

_ __ n SF M P 

° y-~~kL + A (50) 

o cos 

2 

whereas the appropriate probabilistic central safety factor is determined from the 
equation: 
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E&y ) = 


(51) 


s x M 


S cos 


kL 

2 


P 

+ — 


A 


where 


k = 



(52) 


(c) During Use of Element the Concentrated Moment Remains Constant, 
Axial Force Increases 

The deterministic safety factor hsf is determined from the equation 


M 


<7,. = 


n SF P 


S cos 


kL 


(53) 


where 



j n SF P 
' El 


The probability analog of this equation reads: 


E(l y ) = 


M 


^cos 


kyL 

2 


S,P 

+ - ! — 

A 


(54) 


(55) 


where 


k 


V 



(56) 


For example, let M— 2 kN.m, P = 100 kN. The cross-sectional area is annular 
with mean diameter D m =10 cm; thickness = 0.5 cm, L = 3 m. Then the deterministic 
safety factors become, in three settings 
Case (a): ksf = 1 .85 
Case (b): usf= 3.37 
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Case (c): «sf = 2.52 

In probabilistic setting, if E(L } ) = 300 MPa, the same “central 41 safety factor is 
obtained. Yet, straightforward application of the definition E(L)/cr would be incorrect, 
since the load P and the stress a are not interrelated linearly. 

In some cases, a simplified analysis can be performed: Consider the column that 
is simultaneously subjected to the uniform distributed load q and axial compressive load 
P. Then the exact analysis of the differential equation 

EIw"+Pw = M t (58) 


where w = displacement, M, = bending moment due to the transverse load, 


qL q 2 

M, = — x- — x 


(59) 


leads to the maximum bending moment 


A/., = 


, kL 

1 - cos — 

3 _ 2 _ 

k 2 kL 
cos — 

2 


(60) 


where k = -JpTE! . The approximate relationship, as is well known from the applied 
theory of elasticity reads 

M m &± 

max P_ (61) 

Per 

Thus the stresses would read 


a = 


qL 2 ! 8 


P 

■ H 


( 62 ) 


If both q and P increase simultaneously, the safety factor is found from equation 
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n SF qL 2 /8 . n SF P 

y ( n P A 
S 1-^ A 

P. 


The probabilistic “central” safety factor is found from the equation 


rr , s i p 

E(ly) ~ 7 7 p\ ~a 

S 1-^ A 


and not what would appear appropriate at the first glance: 


EVy) 

qL 2 / 8 P 

P^A 


The correct equation (64) leads to a quadratic equation for the central safety factor s\: 


2 SP 2 f SPE(L y ) SP) qL 


-- + SE(L) = 0 


A 8 


Case 5: Yield Stress Has a Weibull Probability Density, Actual 
Stress Is Deterministic 

Consider now the case when the probability distribution function of the yield 
stress reads 


(<r v ) = exp -exp 


% - 


where a z and K are positive parameters. 


Reliability becomes 
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R = Prob(a < E y ) = Prob(I y > a) = 1 -Prob(E y < a) 

= 1-^1, (^) 


( 68 ) 


thus yielding the reliability as follows: 


R = exp 


-exp 


Qy 


V 


/ J 


(69) 


According to Haidar and Mahadevan (2000) (who do not deal with safety factors in the 
present context, but discuss the Weibull distribution) the average and variance of L y are 
directly expressible via a z > and b z > : 


Therefore, 


1 _ 1 n 


a. = ni 1-0 5772*, 


E(l,) = a z +0.5772fc 


Var(Z,)=^b Zy ‘ 


Substitution into Eq. (69) yields 


R - exp 


-exp 


f j6Var(Z y ) 

£'(Z v )- 0.5772— — 

n 


-a 


ij6Var(Z y ) / 7r 


(70) 

(71) 


(72) 

(73) 


(74) 


or 


R = exp 


exp 


' Eq y )-0A5j6Vard y )-<r ' 

0.1Syj6Var(L y ) 


(75) 
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Dividing both the numerator and denominator in Eq. (75) by <x, and recalling the 
definition of central safety factor 5, -E(Z v )/a and variability coefficient of the yield 


stress v £ = ^JVar(L y ) /E(L y ) we get 


R = exp 


^ s, - 0.45v v -1 N 


exp 


v 0.78v z> sj j 


(76) 


If the reliability is fixed, one can find the appropriate safety factor: 


1 + 0.45v,. 


1 1-0.78v e ln(lni) 


(77) 


Since this formula yields safety factors that are less than unity, the use of the Weibull 
distribution appears to be questionable in the case in question. This leads to an all- 
important lesson: Direct randomization of the deterministic problem not always may be 
advisible. 


Conclusion 

In this report we dealt with the case that is reverse to that discussed in report #1. 

As Mischke (1970) notes, “there is disenchantment with the term factor of 
safety.” Shigley (1970) writes: “One of the unfortunate facts of life is that there are 
almost no publications of data on the distribution of stress and strength.” 

We fell that the expert opinions on one hand, and the expert systems on the other, 
along with the use of accumulated data available in engineering firms, will close the gap 
between the present safety factor design and the probabilistic approach to its both 
justification and the rational allocation. 
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Interrelation Between Safety Factor and Reliability: Part 3 


Both Actual Stress and Yield Stress Are Random 

by 

Isaac ElishakofF 

Department of Mechanical Engineering 
Florida Atlantic University 
Boca Raton, FL 33432 


In this part we consider most realistic case when both the yield stress E r and the 
actual stress Z are represented as random variables. The reliability reads: 

R = Prob(E < Z v ) (1) 

We denote by f zl (cr,c x v ) a joint probability density function of 2 and 2 V . Then Eq. (1) 
becomes 


R = jj/n, ia,a y )dada y jf ZZy (a,a y )da y do 


v<v 


( 2 ) 


or, alternatively. 


R = JJ fzz y °y )dcsda y =| 


> 

jf ZZv (<7,a y )da do y 


( 3 ) 


where the integration domain extends over the region in which 2 and E y to be 
independent random variables. We find two formulas, stemming from Eqs (2) and (3), 
respectively: 


R = \\\-Fz y {v)lfx (.a)do (4) 

0 
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(5) 


oo 

R = jF t (<r y )f Cr (<r y )d< T y 


0 

We will use either of Eqs. (4) or (5) based on convenience of computation. 


Case 1: Both Actual Stress and Yield Stress Have Normal 
Probability Density 


Let 


/i(o-) = 


by-Jlir 


exp 


1 f cr-£(X)Y 


fr K) = /= 

' by -Jin 


exp 


(°,-e<x,) Y 


where 

£(X) = mean value of the actual stress 
b z = -JVar(L) — standard deviation of the actual stress 
E(L y ) = mean value of the yield stress 
b Zf = y]Var(L v ) = standard deviation of the yield stress 
We introduce a new random variable 


( 6 ) 

(7) 


M=I V -X (8) 

which is called the safety margin. Since Eq. (8) expresses linearly X and X^, the safety 
margin, as a linear function of the normal variables, is also a normal variable with the 
mean value 

E(M) = £(X V ) -£(X) (9) 
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and standard deviation b m found as follows 


Kt - yjb/ + b z„ 


The reliability is then 


R = Prob(I < S v ) = Prob (M > 0) 


00 1 

f 1 

I py n 


f t-E(M )) 

2 

I , exp 

0 l> M dlx 


l b M J 



Vt 


To perform an integration in Eq. (11), we introduce new variable 

t - E(M) 


Hence 


M 


dt = b M dz 


Also, when t = 0, the lower limit of z equals 

_0-£(M)_ E(l y )-E(L) 


-'M 




Hence, 


1 00 

R - —j= J exp(-z 2 / 2)dz 

<2n e{£,)-e(X) 


Thus, reliability in Eq. (1 1) becomes 


yjb I 


E(Z y yE(Z) 


yjh + %‘ 


R = 1 — p= j" exp(-z 2 / 2)dz 

din L 


= !-<!> 


' £(£„)-£(xf 


Ed 


+ by 


L y J 


( 10 ) 


( 11 ) 


( 12 ) 


(13) 


(14) 


( 15 ) 


(16) 
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This formula can be rewritten in several alternative ways. First of all we note that 


E(l v )-E(l) i 




v m 


( 17 ) 


where v M is the coefficient of variation of the safety margin. Thus 

/? = l-0>[ 


f 1 ' 


V v uJ 


( 18 ) 


We also introduce coefficients of variation of the actual stress and the yield stress, 
respectively, 

E 


Vv =- 


E( I) 

E 

"'y 

E(Z y ) 


( 19 ) 


( 20 ) 


Then, 


4 


b v " by 




2 . V E 2 (h) 


E 2 (L y ) E 2 (I) 


( 21 ) 


= E(Z)Jv ) 2 + v^V 

Hence, the reliability in Eq. (16) becomes 


^ = 1-0 


f 

\ 


5, -1 


/ 2 2 


/v v +V v s. 

V 

V “ -.v 1 J 


( 22 ) 


As is seen the reliability R, the central safety factor si and the variabilities v z and 
are directly interrelated. 
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Case 2: Actual Stress Has an Exponential Density, Yield Stress 
Has a Normal Probability Density 


For the titled case the probability densities of L and ~L y read, respectively, 
/i (p) ~ a exp(-acr), for a > 0 

/i ’ <<7,)= v!& exp 




, for 


00 < (J < oo 

— V 


( 23 ) 


We note that 


£(!) = -, For(2) = 4r 

a a~ 

£(!,) = k, Var(L y ) = b z 1 


We evaluate reliability function as follows 


0 


j /v ( <j)d<7 do 


Now, the inner integral equals: 

a y <T V 

\fzip)dc 7= J aexp(-ao)do 

0 0 

= 1 - exp(-acr v ) 

This results in the following evaluation: 


( 24 ) 

( 25 ) 


( 26 ) 
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R = J y=exp 

l b Y din 

1 00 

= 7= jexp 

by din i 


r a y -E(Z^' 


\ ~y s 

r a y -E( I v ) V 


(1 -e°° y )do y 


by y[ 2 n ‘ 


CO 

r 

1 

r <r, -E(i,)) 


Jexp 

2 

by j 


0 


V z > J 



do. 


e“° y do. 


( 27 ) 


= 1 - 0 ) 




t=fexp -£(Z,) + a»i 2 + 2 E(Z, )b, 2 -a b 

Uy V 27r | dLi)y 


We introduce the following variables 


er v -£(I v ) + a6 r 
r = — : 

by 


Hence 

6 V Jr = Jcr v 

**r ' 

The expression for the reliability reads 


/? = l-0> 

where 

/(/) = expj 
leading to the final expression 


E(L,) 


1 


din 


CO 


E(X y )-aby 2 


V. 




v 


exp|~- — (2f?(£ v )a - a ' ) 


(28) 


(29) 


(30) 


(31) 



£(£ v )’ 


1 o 1 




R = l-0> 

L > J 

-exp 

~-( 2 E(-L y )a-a 2 b z 2 

l-0> 

.v 

\ z y J 



( 32 ) 
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This expression is rewritten as follows with notation: 


v v =• 


E(Z V ) 


( 33 ) 


we get 


exp 




i 

f £(£,,) Ear(I v )y 

-~\2E(L y )a - a~b z ) 

= exp 

2 

l E(L) £ 2 (S) J_ 


exp 


- exp 


- 


Ear(L v )£ 2 (S v ) 
E 2 (E v ) £ 2 (X) 


~VV) 


(34) 


Also, the expression in Eq. (32) becomes 


1-<J> 


E(Z y )-ab z ; 

— 1 — 0 

E(L y ) JVar& y ) 

by 

<JVar(L y ) E(Z) 


= 1-0 

f ] JVarP y ) E(£ y f | 
[ v Ef E(S y ) E(Z) J 

f 1 ^ 

Vy S, 


= 1-0 


(35) 


Thus, all parameters in Eq. (32) are expressed in terms of the coefficients of variation and 
the central safety factor s, 


R = 1-0 


' 1 ' 




-exp! 


1 / 2 2 \ 
--( s i ) 



u 

** 

i 

i 


1 - d> 



i j 



(36) 


Consider an example. Let the yield stress has a normal probability density with mean 
yield stress equal £(£ v ) = 100 MPa. The variance equals Var(Z ) = lOO(MPa) 2 , or 


standard deviation equals b z =10 MPa, leading to the coefficient of variation to be 
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(37) 


10 


Vy = 


E(Z V ) 100 


= 0.1 


The central safety factor is set at 2, i.e. 

EV v ) ~ 

s, = — = 2 

E(l) 

leading to the value of the mean stress to be equal 

E(L) = —E(L V ) = — x 100 = 50 MPa 
2 2 


(38) 


(39) 


i.e. 

a = \/E(L) = 1/50 (MPa)' 1 (40) 

Calculations in accordance to the formula (36) yield the reliability R = 0.86194. 


Case 3: Actual Stress Has a Normal Probability Density, 
Strength Has an Exponential Probability Density 


In this case we get 


00 co 

R=\M<T) J h s {° y )do y do 


(41) 


In new circumstances the probability densities read: 


/i(ct) = 


^2nVar(T) 


exp 


4K)= 


EV v ) 


exp 


E& v ) 


o - £(S) 
yJVar(Z) 


|> <^v > 0 


(42) 


The reliability becomes: 
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M- 


f 1 



2“ 

exp 

/ \ 

a 

’^InVari!) 


{ 2 Var(L) j 




d<7 


(43) 


The find expression is as follows: 

R = -X-, , . , -,2 E(S) Var(D 


E(l) ^ 


r if 

yjVar(L) j 

+ exp 



y 

f 


1-<I> 

)_ 

V 


( £(£)-Far(I)/£ 2 (S v ) 
^Var(Z) 


E&y) E 2 Vy) 

In terms of coefficients of variation and the central safety factor this expression becomes. 


(44) 


R = t D 


0 


i 

f ~ 2 V 

2 Vy 


r _J__ VA 

— 

+ exp 

— 

2 

1-0 




2 

1 5 i s \ J . 


l v s /_ 


(44) 


Case 4: Both Actual Stress and Yield Stress Have Log-Normal 
Probability Densities 

Let the actual stress have the following density 




aby^lrr 


exp 


(In cr -a z ) 2 
2b 2 


, cr > 0 


where parameters a z and b z are related to the mean value E(L) as follows: 

E(Z) = exp (a z + '-b z 2 ) 

The variance of the stress equals 

Var(Z) = e\p(2a z + £. 2 )[exp(Z> £ 2 ) -1] 

The probability density of the yield stress reads: 


/i.K) = 


1 


cr v bv >f2n 


exp 


(In <7 y -a £> ) 2 
2b7 2 


, <r y >0 


(45) 


(46) 


(47) 


(48) 


The parameters a z and b z are related in the following way with the mean yield stress: 
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E(I V ) = exp(a^ +\b t> 2 ) 


(49) 


Variance Var(L y ) is expressed as 

Var(Z v ) = exp(2a z +b z ')[exp(6 z ')-!] 


The reliability reads 


R = Pr ob(I < E ) 


Yet, it is easier to express reliability as follows. 


fl = Prob 


A y 


Introducing a new variable Z, 




we get 


R = Prob(Z < l) 


which can be rewritten as follows: 

/? = Pr ob(ln Z < 0) 


We note that 


InZ = InE- ln£. 


But In I has a normal probability density with 

£(lnE) = a s 

Var( lnl) = V 

Likewise In Z y has a normal probability density with 


(50) 

(51) 

(52) 

(53) 

(54) 

(55) 

(56) 

(57) 

(58) 
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£(lnE v ) = a Zt 
Var( In S v ) = b z 

Hence, the difference In Z = In L - In Z v has a normal probability density with 


( 59 ) 


E(\nZ) = a z -a z 


Par(lnZ) = + b z 


(60) 


We also conclude that Z is log-normal random variable with 

E(Z) = exp[a s - a £ + y (& £ 2 + b z " )] 

Var(Z) = exp[2(a £ -a £ ) + 6 S 2 + b z 2 ][exp(A s " + b z ~ ) - 1] 


The reliability reads 


R = m 


a z -a t 


V 


^ 4" 


(61) 

(62) 


(63) 


The central safety factor reads 


E(Z. y ) exp(a £ + b z 12) 


E(L) exp (a s +b 1 ~/2) 


(64) 


Consider an example. Let 

E(L) = 60,000 kPa 
yjVar(Z) = 20,000 kPa 
E(Z V ) = 100,000 kPa 

yfVar(L^) = 10,000 kPa 

We know that the central safety factor equals 

£(£ v ) 100,000 


=' 


£(£) 60,000 


= 1.67 


(65) 


( 66 ) 
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The reliability in this case turns out to be equal R = 0.9495. 


Case 5: The Characteristic Safety Factor And the Design Safety 
Factor 


The characteristic safety factor reads 

y 

" v.0.05 

7 ~ — — 

y 

" 0.95 


(67) 


where Z v005 is 0.05 fractile of the probability distribution of yield stress, Z 095 is 0.95 

fractile of the probability distribution of stress. 

When both random variables I and Z>, are normal, then, according to Leporati [2] 

£(Z V )(1 -1.645v r ) 1- 1.645V,. 

/ = — = 5 , — ( 68 ) 

E(Z)(1 + 1.645vv) 1 + 1.645v s 

If, for example, both coefficients of variation are set at 0.05 then characteristic safety 
factor equals 


0.91775 
1.08225 i] 


0.848.9, 


(69) 


For the coefficients of variation set at 0. 1 we get 

y = 0.316^ 

The design safety factor, according to Leporati [2] equals 

v 

* ^ v,0 005 

7 ~ — 

y 

"0.95 


(70) 


(71) 


where Z 0005 is the 0.005 fractile of the yield stress. It equals 

l-2.576v^ 

7 ~ Sx 1 + 1 ,645v e 


(72) 
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Case 6: Asymptotic Analysis 


Eq. (16) can be put in the following form: 

P f = 4*0) 


( 73 ) 


where flis referred to as a reliability index. 

For fi> 5, the probability of failure can be written via an asymptotic formula: 

1 p 2 -\ f p 2 ) 

t) 

The reliability index itself is represented as follows, via Eq. (16) 


(74) 


£(£',)-£(!) 
JVar(I v ) + Var(l) 

Dividing both numerator and denominator by E(L) we get 

p=- r J±= L= 

^ / -> 2 2 

JVv‘ + 5, V z> 

From Eq. (76) we can find S\ by solving the quadratic 

p 2 v z + PW 2 = (•*, -l) 2 
s 2 (P 2 Vy 2 -1) + 2 s ] + p 2 v z 2 -1 = 0 


(75) 


(76) 


(77) 

(78) 


We get 


1+ V 

h-iPW - \)(P 2 V Z 2 - 1) 


l~0\ 2 

1 + \ 

lp\ 2 +fiW -0 4v z,\ 2 


1 -p 2 Vy 2 


(79) 
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As is easily seen, not for all coefficients of variability v z and v z> one can find the 

central factor of safety such that the demanded reliability index would be achieved. For 
example for 


(fi 2 v z 2 -l)(j3 2 v z 2 -1)>1 (80) 

5] gets complex values 

Solving the inequality (80) we get 

P 2 (P\ 2 Vz , 2 -v s / ~V)>0 (81) 

This implies, that the following inequality must be met 


P< 



(82) 


If the inequality (82) is violated, then the required reliability cannot be achieved by any 
factor of safety. 

From Eq. (76) we observe that when Si tends to infinity, the reliability index tends 
to 


Vv 


(83) 


We can differentiate Eq (76) with respect to Si: 


V £ +V X„ 


dl = 

ds } (v t 2 + s, 2 v £ 2 ) 


2 \ 3/2 


>0 


(84) 


This implies that when increasing 5 , from unity (when fi= 0) to infinity, ft varies from 
zero to the value 1 / v r . 


If the variation of the stress is zero v z = 0 Eq. (79) yields 
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1 


( 85 ) 


1 -J3v z 

When we have a zero variability of the yield stress = 0 , we get 

s i ~ 1 + A’v 


( 86 ) 


Case 7: Actual Stress and Yield Stress Are Correlated 


“In most cases the correlation between the actual stress and the yield stresses is 
absent. Yet, when such a correlation exists, its expressing is ambiguous, and it is difficult 
to express it numerically, “as Rzhanitzin (1981) notes. 

Positive correlation between I and 2y take place when the stronger elements take 
more load. Partially this takes place for statically indeterminate systems, in which greater 
strength is associated with greater stiffness, and hence with more loads. Safety margin 

M = 2 V -Z (87) 

has a variance 


Var(M) = Var(L v ) - 2 Cov(Z 2) + kar(S) 


( 88 ) 


where Cor(2 ,2) is the covariance between the actual stress and yield stress. 

Then instead of Eq. (75) we get 

E(LJ-E(L) 


P = 


JVar(L y ) - 2Cov(I v , 2) + Var(Z) 


(89) 


Hence, in terms of central safety factor, we have 




s x -1 


f 


£ 2 -2*1 V EI 2 +5 l 2y £ 2 


(90) 
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where v„ is the correlation coefficient 


I Cov(Z v ,Z) 

V £(X y )E(Z) 


(91) 


One can express s\ via p Eq. (90) becomes 




(1-Vj. 'P )$," -2(1 - P L \\v )s t +(l~P~Vy-) = 0 


(92) 


yielding 


J i = 


1 P v z y + P ( v ij, + v s ) P ( v i V i v v n: J , ) 


(93) 


When Vq. vanishes Eq. (93) reduces to the case of uncorralated actual stress and yield 
stress. 


Case 8: Both Actual Stress and Yield Stress Follow the Pearson 
Probability Densities 


The random variable is said to have a Pearson probability density, if it has a form 
f x (x) = Ax a e~ bx for x>0 (94) 

where 


Ha + 1) 

The mean value E(X) reads 

E(X) = ^± 

whereas the variance equals 


(95) 


( 96 ) 
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Var{X) = ^± 
0 


( 97 ) 


The coefficient of variation equals 


Vv = 


(98) 


- v VoTT 

It is interesting that the coefficient of variation does not depend upon b. It is seen that for 
small values of a we get very high variability, whereas for small variability, of the order 
of 0. 1 , greater values of a are needed (a * 1 00). 

Let the actual stress have the Pearson density 


/ s ( o -) = A °‘“" 1 ex P 


acr 


V E(T)J 


The yield stress also has the Pearson density 


where 


h (<t v ) = A 2 a/~‘ exp 


P°y 

£(£>.) 


(99) 


( 100 ) 


where 


A _ - A - £1 

1 r(a)£ a (S)’ 2 T(P)E p {Y) 


a = -K, /? = — 

v y " v y 


(101) 


( 102 ) 


The probability of failure becomes 

p B s (0,a) 
f B(P,a) 


= J s (a,P) 


(103) 


where 
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B ^ a)= n a w m 

T{a + fi) 


( 104 ) 


is the Euler function of the first kind, or beta-function, whereas 

s 

B s {p,a) = \t p \\-ty- l dt (105) 

0 

is the Euler function of the second kind, or incomplete beta-function. The quantity 8 in 
Eq. (105) is defined as follows: 


/?£(£) 

/JE(I) + aE(l y ) 


1 + 



(106) 


The evaluation for the function can be done by the numerical evaluation of 

integrals in Eqs. (105) and (106). 


Conclusion: Reliability & Safety Factor Can Peacefully Coexist 


At this junction we ask the most important question: Are the safety factor and the 
reliability concepts contradictory or they can coexist peacefully? The Fig. 1-12 depict the 
dependence of the reliability in Eq. (22) vs the central safety factor sj. As is seen, for 
various variabilities of the stress and the yield stress one can assign both reliability and 
the central safety factor. Indeed, Fig. 1 shows that of the coefficient of variability of the 
stress ye= 0.04 and the designer wants the safety factor to be set, say at 1.3, the demand 
that the reliability is above 0.92 results in the choice of materials with yv y <0. 15. We 
immediately observe that both the central safety factor and the reliability requirements 
can be combined. Likewise, for yv=0.05 (Fig. 2), the central safety factor s=1.2 is 
associated with reliabilites grater than or equal to 0.93. Thus, one can also impose the 
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reliability constraint. Analogous features are characteristic to figures 3-1 2. This leads to 
the conclusion that these 2 concepts can coexist peacefully. 

But this coexistence cannot be done without some adjustments Reliability 
concept provides the rigorous values of the factors, for otherwise, i.e. without the general 
context of reliability, the safety factor will remain as the factor of experience but still the 
factor of theoretical ignorance. Adopting reliability as the main concept, the allocated 
values of safety factor will naturally follow. In a sense probabilistic methods do not 
constitute a “revolution” but rather a natural “evolution”. 
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Fig 3: Reliability, R, versus Central Safety Factor, s 
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Fig 5: Reliability, R, versus Central Safety Factor, s 
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